Evaluation of ranks of real space and particle entanglement spectra for large systems 
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We devise a way to calculate the dimensions of symmetry sectors appearing in the Particle Entan- 
glement Spectrum (PES) and Real Space Entanglement Spectrum (RSES) of multi-particle systems 
from their real space wave functions. We first note that these ranks in the entanglement spectra 
equal the dimensions of spaces of wave functions with a number of particles fixed. This also yields 
equality of the multiplicities in the PES and the RSES. Our technique allows numerical calculations 
for much larger systems than were previously feasible. For somewhat smaller systems, we can find 
approximate entanglement energies as well as multiplicities. We illustrate the method with results 
on the RSES and PES multiplicities for integer quantum Hall states, Laughlin and Jain composite 
fermion states and for the Moore- Read state at filling v = for system sizes up to 70 particles. 



PACS numbers: 71.10.Pm, 73.43.-f, 05.30. Pr, 11.25.Hf. 

In recent years, entanglement has been embraced as 
an important tool in the study of many particle systems. 
After splitting the system into subsystems A and B the 
entanglement between these can be quantified using the 
entanglement entropy Sa = — tr(yOA log(pA)), where pA is 
the density matrix of subsystem A. One splitting is ob- 
tained by dividing the space in which the system resides 
into disjoint parts. For topological phases of matter, one 
finds that Sa then increases in proportion to the size of 
the boundary between these, up to a constant topolog- 
ically invariant correction, the topological entanglement 
entropy 7 PQ [2] . In an attempt to calculate 7 for frac- 
tional quantum Hall (QH) states, the authors of Refs. 3 
and introduced a splitting of Hall systems in angular 
momentum space, the orbital cut. Since the orbitals of 
the Landau problem localize the particle near rings or 
lines whose location relates to the angular momentum, 
one may view the orbital cut as a cut in real space. How- 
ever, as spatial regions of high probability density can 
overlap strongly for different orbitals, the orbital cut is 
actually very 'fuzzy' in real space. A different cut which 
simply divides the particles in two groups, irrespective 
of their positions, was also considered in Refs. [SHSl We 
will see that, perhaps counterintuitively, this cut relates 
more directly to a sharp spatial cut. In Ref. 121 Li and 
Haldane introduced the entanglement spectrum as a gen- 
eralization of the entanglement entropy. They consider 
all eigenvalues of pA, and the corresponding states in the 
Schmidt decomposition. Li and Haldane conjectured a 
correspondence between the states in the orbital entan- 
glement spectrum (OES) of the ground state of a Hall 
system and the edge excitations of one of the subsystems. 
This implies that one may make study of the edge exci- 
tations of many gapped systems even if only the ground 
state of the system can be accessed directly. A simi- 
lar conjecture states that the particle entanglement 



spectrum (PES) should reproduce the spectrum of the 
smaller subsystem at the same total flux as the full sys- 
tem. These correspondences between energy and entan- 
glement spectra are necessarily qualitative, since relevant 
data like edge potentials are not represented in the bulk 
ground state. Also, in model systems, some classes of 
excitations may not be reflected in the ground state - 
one may think of fully spin-polarized systems which have 
excitations involving spin flips. Nevertheless, generic ar- 
guments in favor of the OES-edge correspondence were 
presented in Refs. |5] and [HI in cases where the edge has a 
conformal field theory (GET) description. In Ref. 10. it 
was also shown that, for a large class of Hall states, the 
PES contains the same information about the edge as the 
OES in the thermodynamic limit. Moreover, the corre- 
spondence holds for systems beyond QH states, including 
topological insulators and superconductors [TTl [T^] . 

Here, we calculate the angular momentum multiplici- 
ties of the real space entanglement spectrum (RSES), i.e. 
the ES which results from a sharp cut in real space (as in 
the definition of the topological entanglement entropy). 
These multiplicities are the ranks of the angular momen- 
tum blocks of the reduced density matrix. For QH states, 
the RSES and RSES ranks have so far been calculated 
only for the non-interacting states at some integer fill- 
ings [T31 [T3] . The method of calculation that we develop 
works just as well for the PES ranks and we will soon 
see that the RSES and PES ranks are in fact equal. This 
means in particular, using the results of Ref. I10[ that, in 
many cases, the RSES ranks contain the same informa- 
tion about the thermodynamic limit of the edge as the 
OES ranks, though finite size effects are different. In gen- 
eral, we expect the RSES to give a better description of 
the edge than the OES. For example, for the integer Hall 
states, the RSES ranks reproduces the edge multiplicities 
faithfully, while the OES is trivial. 
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Consider a system of N particles in its ground state 
tfj. We start with a general statement of the conjectured 
correspondence between the entanglement spectrum and 
energy spectrum of such a system. Assume that the sys- 
tem has a symmetry J = X^iLi Ji ^^'^ ^^^''^ V' is an eigen- 
state of J. If the splitting of the system preserves J, 
then the reduced density matrices of the subsystems will 
be block diagonal, with the blocks labeled by the eigen- 
values of J A (or Jb = J — Ja)- It is natural to expect 
that the number of (low lying) states in the entangle- 
ment spectrum obtained from the ground state '0, with 
a given eigenvalue of J a, will be the same as the number 
of (low energy) states of subsystem A which have that 
same eigenvalue of J a- One may further hope that the 
entanglement energies and the real energies of the low- 
lying states become proportional in the thermodynamic 
limit and that the Schmidt states approach the low lying 
states in the energy spectrum. Here, we focus mainly on 
the ranks of the blocks of pA, but we do present prelimi- 
nary results on entanglement energies (see Fig. |2| . 

Now take i^^zi, .., z^) to be a QH wave function (with 
Zi the complex coordinates of the particles) and take J 
to be the angular momentum component L^. The sub- 
systems A and B contain k and N — k particles. The 
matrix elements of the blocks of pA are then 

p^^'\zi...Zk,z[...z',^) = J nf=fe+i dzj 

J d</.e-2'^*^^'^V(2ie^''^'^ . . . ^fee^'^*'^, Zk+i . . . zjv) 
/ d(j)' e^^^LW^*{^z[e^^"!'' . . . z^e^'^*'^' , Zk+i ...zn). 

(1) 

Here the role of the (j) and 0' integrations is to enforce 
that the angular momentum in subsystem A equals L^, 
therefore selecting a single block of pA- Defining Za — 
(zi,..,Zfc), Zb = (zfc+i,..,2;jv) and 

$(L^,Z^,Zb) - 

j #e-2^'^^*V^(zie2"^..,Zfee2'^^^Zfc+i,..,zw), (2) 
we can rewrite ([T]) as follows: 

p^j:^\ZA\Z'j,) =jdZB<^{LtZA,ZBW{LtZ'A,ZB)m 

If the positions in Za can take any value and the integrals 
over Zb are over all space, then the above expressions 
are for the reduced density matrix with the particle cut. 
However, we can restrict the Za positions to lie in a 
subset of space (e.g. a disk) and take the integrals over 
Zb to be over the complement of that subset. The same 
expressions then represent density matrix blocks for the 
real space cut. We can even allow the regions containing 
the Za and Zb positions to overlap, giving cuts which 
interpolate between particle cut and real space cut. 

The right hand side of ([S]) is the Gram matrix for a set 
P of wave functions for the particles in subsystem B: 

^ = {^L-,zAZB) = HLtZA,ZB)}. (4) 



Therefore, the rank of pA is equal to the number of lin- 
early independent wave functions in P. The wave func- 
tions in P are precisely the = — components 
of the wave functions for a system of Nb electrons in 
a background of Na fixed electrons at positions Za, or 
if you will, for Nb electrons with Na 'electron holes' 
at positions Za- While the positions Za are continuous, 
the number of linearly independent wave functions of this 
kind is finite as long as the Hilbert space for Nb particles 
contains only a finite dimensional subspace at the given 
angular momentum. This is true for wave functions built 
from a finite number of orbitals, e.g. from orbitals in a 
finite number of Landau levels. Much stronger bounds on 
the ranks hold when ij) satisfies vanishing properties when 
some number of electron coordinates coincide (e.g. if ^ is 
the bosonic Moore- Read (MR) Pfafhan [TS], it vanishes 
whenever 3 particles coincide). Since all wave functions 
in P satisfy the same vanishing properties, the ranks are 
bounded by the numbers of independent states with the 
given vanishing properties at the given magnetic flux. 

We now argue that the ranks of the blocks of pA are the 
same for the particle cut and real space cut. In both cases 
the ranks equal the numbers of independent wave func- 
tions in P. The functions in P actually remain the same, 
only the domains on which they are defined are different; 
in one case the positions Za and Zb can be anywhere in 
space and in the other case, they are restricted to two dis- 
joint regions. However, up to geometrical factors, these 
functions are polynomials. If two polynomials are equal 
on any open domain, then by analytic continuation, they 
must be equal everywhere. Therefore the number of inde- 
pendent states is the same for particle cut and real space 
cut. We should stress that despite the equality of ranks, 
the actual eigenvalues of pA will be very different for the 
two cuts and there is physical information contained in 
these. For example for the RSES we expect lower entan- 
glement energies for states associated with the edge. In 
fact, we could choose the cut so that subsystems A and 
B share multiple circle boundaries and in that case we 
would expect that each edge will give rise to a branch of 
low entanglement energy states in the RSES. 

Given trial wave functions for quasiholes, we can at- 
tempt a direct proof of the correspondence between the 
number of independent wave functions for quasiholes and 
the entanglement spectrum rank by expanding the wave 
functions ^l^.Za an appropriate polynomial basis. If 
the same count is obtained for the ^j^a zj^ as for quasi- 
hole wave functions at the same flux, then the correspon- 
dence is proved. Note that the quasihole wave functions 
which should be compared with the ^l'^.Za ^^^^ have 
-^Na quasiholes, where e* is the charge of a single quasi- 
hole. We will not pursue this here. Instead, we develop a 
numerical technique for the calculation of the RSES and 
PES ranks, which can be applied for any ■(/;. 

Our basic strategy is as follows. Even though the set 
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of index configurations Za^'Za' is infinite, the rank of 

Pa is finite. Therefore we can find it by evaluating the 
rank of a suitably large submatrix which we obtain by 
choosing a particular set of values of Za and Za'- We 
will use a square d x d submatrix with the same set of 
Za and Za' , which we call Z\ with i £ 1, . . . , d. The 
rank of this matrix is the same as the number of linearly 
independent wave functions in the set 

P = {^,{ZB)^HLtZ\,ZB) i = l,..,d}. (5) 

This number is equal to the rank of the matrix 

= l>,(Z^) z = 1, . . . , d, j = 1, . . . , d' > d (6) 

where Z^ = (z^+i, ■.ZN)j {j = 1, ..,d') is a set of d' > d 
different (N — fc)-tuples of coordinates in the B subsys- 
tem. The rank of this matrix can be obtained (e.g. by 
singular value decomposition) and equals the number of 
nonzero eigenvalues of M. In numerical calculations, 
all eigenvalues of M^M will be nonzero, but a clear jump 
between large and small (nearly zero) eigenvalues is ob- 
served and hence the true rank can be read off. 

The scheme just described will not work without a ju- 
dicious choice of the index configurations Z\ and Z;^. 
Choosing these at random will lead to very small and 
greatly varying values of the matrix elements . This 
induces numerical error which makes it difficult to iden- 
tify a clear cut in the spectrum of eigenvalues of iW^M. 
A better set of index configurations can be obtained by 
Monte Carlo sampling ip, but harvesting only configura- 
tions whose Za satisfies a suitable constraint, which helps 
to select index configurations which have good overlaps 
with the desired angular momentum sector. The con- 
straint is introduced by means of the relationship be- 
tween the coordinates and the angular momentum in the 
Landau levels. For instance on the disk the single par- 
ticle state with angular momentum I has most of the 
probability density concentrated near a ring of radius 
Therefore, we can assign to each Metropolis con- 
figuration Za = {zi, .., Zk) a A;-tuple of angular momenta 
Ia = {Ia: --i^a) by taking l\ equal to the integer nearest 
to I Zip. We then harvest only those Za which satisfy 
J2i=i — ^ index configurations. 

When two configurations Za and Z'j^ obtained in this 
way share the same I a , we can discard one of them with- 
out loss of rank, because the corresponding rows of the 
matrix M are almost exactly proportional. Moreover, to 

obtain the full rank of ° it is not necessary to gener- 
ate index configurations corresponding to each admissible 
/c-tuple I A- This is because the set of states labeled by 
the Za is overcomplete and states with different values 
of Ia have nonzero overlaps. To obtain the full rank r 

of one of the , we typically need to take a num- 
ber of indices d which is only a few times larger than r. 
If we want to focus on the low lying part of the RSES, 
a further simplification takes place. Now the electron 



coordinates in the Z\ must also lie inside the spatial re- 
gion associated with subsystem A. E.g. the k particles 
in subsystem A may be located in a disk of radius ta, so 
that \zi\ < Ta for i E {!,.., A;}. In this case the associ- 
ated angular momenta l\ must also have their probabil- 
ity density concentrated on rings located within this disk 
and the number of Ia which satisfy the constraint that 
X]i=i 'a = -^f is much smaller than if zi, . . . ,Zk could 
be located anywhere. This selection of indices allows for 
a very efficient evaluation of the RSES ranks in the an- 
gular momentum sectors associated with the edge, but it 
does make it more difficult to find the multiplicities for 
other values of Lf associated with bulk excitations. To 
obtain those it is better to select indices Z\ with electron 
coordinates which vary through all space. 

Using this method, we can manage systems with up 
to iV = 100 and Na = 50 particles on a standard lap- 
top. Note that for systems where the correspondences 
between the ES and the bulk and boundary excitations 
hold, this gives us information on the excitations for sys- 
tems of up to 50 particles, which is considerably larger 
than the system sizes accessible by exact diagonalization 
(typically no more than 20 particles). With more effort 
and resources, much larger systems should be accessible. 
It should also be clear the method can be used on the 
plane, the cylinder and the torus. 

We now present a summary of results of sample PES 
and RSES rank calculations. These were done for integer 
filling fractions on the cylinder and for the Moore-Read 
Pfafhan (TFj , Laughlin and Jain states both on the cylin- 
der and in the spherical geometry [16]. We consider the 
case Na — N/2 and for the RSES we define the subsys- 
tem A as half of the sphere or cylinder. 

For the Laughlin states at = 1/2 and v = 1/3 we 
calculated the RSES ranks for systems of 70 particles. 
The multiplicities of the low-energy edge excitations of 
the Laughlin state, in the thermodynamic limit, are pre- 
dicted by the edge conformal field theory (CFT), which 
is a chiral Luttinger liquid. They are 



edge counting 



0123456 7 8 n 
1 1 2 3 5 7 11 15 22 p{n) 



(7) 



Here AL^ is the relative angular momentum between the 
quasihole excitations and the Laughlin groundstate and 
p{n) is the number of partitions of n. We have computed 
the RSES ranks up to Lf = 20 and find that for both 
h' = 1/2 and u = 1/3, they match the edge counting 
above. Fig. [l] shows the spectra of the matrices M 
for the Laughlin state at v = 1/2 at different values of 
, in the way in which entanglement spectra are usually 
presented {Lf is given relative to the state with the low- 
est Lf, analogously to AL^ above). A clear gap is visible 
for each value of Lf and the counting of large eigenvalues 
matches the partition numbers. Of course this plot is not 
a plot of the true RSES. Repeats of the RSES rank calcu- 
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FIG. 1. Plot of the spectra of the matrices M^M used in the 
calculation of the RSES ranks for a u = ^ Laughlin state at 
N — 70 and Na = 35 on a sphere. The negative logarithms 
of the eigenvalues for each block of pA are plotted (the 
largest eigenvalue is normalized to 1 at each L^). A clear 
gap is visible for each value of and the counting of large 
eigenvalues (at low e) matches the partition numbers. 



lation with different initialization of the random number 
generator will yield different, but similar looking plots. 

For the bosonic and fermionic MR states, we also com- 
puted the RSES at N = 70 and if < 20, obtaining the 
expected CFT countings [17]. We also checked that, by 
changing or by including localized quasiholes, one 
may obtain the CFT countings for different topological 
sectors, in analogy to the results for the OESp, IT5]. 

We have computed the PES ranks beyond the universal 
CFT numbers, for v = 1, i' = 1, and v = \ and for the 
fermionic and bosonic MR states, for many values of 
at N = 'aQ and for all in smaller systems, obtaining 
the expected finite size countings [HI [20] . 

At integer filling g N the OES ranks are trivial but 
the RSES ranks are the full ranks of the angular momen- 
tum sectors of the Hilbert spaces for Na particles in v 
Landau levels. For v = 1, the ranks are just the partition 
numbers (j?]) in the thermodynamic limit, in agreement 
with Ref. [IB] , which treated the v = 1 system analyt- 
ically in a second quantized formulation. For general 
integer filling v = p, the ranks depend on the residue of 
Na modulo p. E.g. for = 2, a generating function for 
the thermodynamic countings is 

VmeZ / \*:>0 ^ / 

where s = A^/i(mod 2). Similar formulas may be obtained 
for higher integer v. In fact, it can be shown that for 
Na = (mod v) the ranks always equal the numbers 
of ^-colored generalized Frobenius partitions [21]. We 
checked that our method reproduces these numbers for 
V < 3, for systems of around 30 particles and necessarily 
for modest angular momenta, as the ranks grow quickly. 




FIG. 2. Plot of the spectra of the matrices Af^M used in 
the calculation of the RSES ranks for a v = ^ Jain state 
on a sphere at iV = 20 and Na = 10 (left) and at iV = 22 
and Na = 11 (right). Negative values of e occur because 
we used an unnormalized Jain state (normalization shifts the 
vertical axis by a constant). The counting 1, 2, 5, 10, . . . of the 
branches is consistent with a pair of noninteracting Luttinger 
liquids. Differences in total counting for even vs. odd Na arise 
analogously to those for u = 2 (cf. Eq. (|8|). 

Finally, we apply our method to the Jain state at 
J/ = |. Jain states [H] have so far produced serious 
challenges in ES calculations because they lacked a clear 
'entanglement gap' at the accessible sizes. While the nu- 
merical method described here calculates RSES ranks, 
we can expect the low lying parts of the spectra of the 
matrices M^AI to converge to yield the low lying parts of 
the true RSES, as the size of the matrices M is increased. 
Fig. [2]shows the spectra of the matrices M^M at low 
for the = I CF state at TV = 20 and TV = 22, ob- 
tained with matrices M which are large enough to allow 
the pattern of entanglement energies to emerge. Note 
that so far, no ES for this state has been published for 
N > 10, due to the difficulty of obtaining the wave func- 
tion in angular momentum space at large N. At low Lf 
values, we observe several branches of edge states in the 
RSES. Each branch has ranks 1, 2, 5, 10, . . ., consistent 
with the counting for a pair of non-interacting Luttinger 
liquid edges. The total rank at each equals the di- 
mension of the spaces of CF states with Na CFs in 2 
CF Landau levels at the given value of Lf. The ranks 
are in general smaller than those for i/ = 2 because some 
CF states disappear in the lowest LL projection. Note 
that the CF LLs we used for the A system are the same 
as those for the full system. This is natural for the low 
edge states of the RSES. For a description of the 
full PES, one should consider CFs at a higher effective 
flux, see Ref. [7] for details. We further checked that the 
space of eigenvectors of M'^ M at each Lf has very close 
to unit overlap with the space spanned by the vectors 
Vi = (f>f^{Zj^), where the (pf^ are the CF trial wave 
functions for system A and the are the configurations 
of the particles in subsystem A that were used in the 
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construction of M. It appears that the low lying RSES 
of this CF state is excellently described by the CF trial 
wave functions for the edge excitations of the A system. 
In future work, we intend to expand our results on CF 
states and on the calculation of entanglement energies. 
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Note added: As this work was being completed, we 
became aware of research studying the RSES by A. Ster- 
dyniak, A. Chandran, N. Regnault, B. A. Bernevig and 
P. Bonderson [53], and by J. Dubail, N. Read and 
E. H. Rezayi [H]. 
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